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ABSTRACT 

We present a chemical abundance analysis of the secondary star in the black hole binary 
V404 Cygni, using Keck I/HIRES spectra. We adopt a x^-minimization procedure to derive 
the stellar parameters, taking into account any possible veiling from the accretion disk. With 
these parameters we determine the atmospheric abundances of O, Na, Mg, Al, Si, Ca, Ti, Fe, 
and Ni. The abundances of Al, Si, and Ti appear to be slightly enhanced when comparing 
with average values in thin-disk solar-type stars. The O abundance, derived from optical lines, is 
particularly enhanced in the atmosphere of the secondary star in V404 Cygni. This, together with 
the peculiar velocity of this system as compared with the Galactic velocity dispersion of thin-disk 
stars, suggests that the black hole formed in a supernova or hypernova explosion. We explore 
different supernova/hypernova models having various geometries to study possible contamination 
of nucleosynthetic products in the chemical abundance pattern of the secondary star. We find 
reasonable agreement between the observed abundances and the model predictions. However, 
the O abundance seems to be too high regardless of the choice of explosion energy or mass cut, 
when trying to fit other element abundances. Moreover, Mg appears to be underabundant for 
all explosion models, which produces Mg abundances roughly 2 times higher than the observed 
value. 



Subject headings: black holes: physics — stars: abundances — stars: evolution 
V404 Cygni (=GS 2023-1-338) — supernovae: general — X-rays: binaries 



stars: individual 



1. Introduction 

The low-mass X-ray binary (LMXB) V404 Cygni 
is one of the most compelling cases for a black 
hole accreting mass from a low-mass compan- 
ion (secondary star). The mass function is 
among the highest values in such binary systems 



syi 

(/(M) = 6.08 ± 0.06 M0; ICasares et al]ll992 : 



Casares fc CharleslHool . placing the mass of the 
compact object comfortably above the upper limit 
of a rapidly rotating neutron star. The mass ra- 
tio of the system, q = Mz/Mbh - 0.060ig:gS35, 
was derived from the measurement of the ro- 
tation velocity of the secondary star, v sin i — 
39.1 ± 1.2 kms^i (jCasares fc Char les 1994). 
Shahbaz et al. ( 19941 ) modelled the ellipsoidal 
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variations using a iC-band infrared (IR) light curve 
of V404 Cygni and derived an orbital inclination 
of i = 56° ±4°, and consequently the implied mass 
of the black ho l e is M rh = 12 ± 2 Mq. Later, 



Pavlenko et al.l (|1996[ ) studied the i?-band light 
curve to determine an inclination of i = 56° ± 2° . 



More recently, iKhargharia et al.l (|2010t ) derived a 
K3 III spectral type for the secondary star from 
near-IR (NIR) spec troscopy of V404 Cygni, and 
Sanwal et al.l ( 1996 ) remodeled the iJ-band light 
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curve; they refined the orbital inclination to be 
which yields a black hole mass of 



I = 



+0.2 



^^o.t M0. 



The system is placed near the G alactic plane, 



with a Galactic latitude of b ; 



2.1°. Miller- Jones et al 



(|20Q9br ) recently derived the distance with very 
long baseline interferometry (VLBI) observations, 
d = 2.39 lb 0.14 kpc. They revisited the analysis 
of Miller- Jones et aP ( 2009al ) to determine a more 
accurate pecuhar velocity of 39.9 ± 5.5 kms^^. 
This value is significa ntly lower their previous 
one, Upcc ~ 64 kms~^ (IMiller- Jones et al. 2009a ) 
assum ing a distance of 4 kpc (IJonker fc Nelemans 
20041 ). and it can be achieved via a iBlaauwl (119611 ) 
velocity kick; thus, no asymmetric supernova kick 
is required. However, Miller- Jones et al.] ( 2009bh 
argued that the component of the peculiar velocity 
in the Galactic plane, 39.6 km s~^, is still too large 
to come from the intrinsic velocity disp ersion of 
the G alactic plane, which is 18.9 km s~"'^( Mignardl 
2fl0fl[ ) for the likely F0-F5 progenitor of the donor 
star. Thus, this peculiar velocity should come 
from a natal kick from the ejection of material in 
a supernova (SN) or a more energetic hypernova 
(HN) event, but its magnitude does not require 
an additional asymmetric kick. 

The chemical abundances of secondary stars in 
black hole and neutron star X-ray binaries have 
been studied for several systems: Nova Scorpii 



1994 ( Israelian et al ■Il999l:l iGonzalez Hernandez et al 



2008aD . A0620-00 
20041) ■ Gentaurus X-4 



Gonzalez Hernandez et al 



Gonzalez Hernandez et al 



2005 ) , XTE J1118-b480 (IGonzalez Hernandez et al 



2006[lGonzalez Hernandez et al.l2008bl). and V464 1 



Sagittarii (jOrosz et al.ll200lHSadakane et al.ll2006l ) . 

The metallicities of these binary systems are all 
close to solar independent of their location with 
respect to the Galactic plane. In addition, the 
above authors have taken into account different 
scenarios of SN/HN ejecta pollution on the pho- 
tospheric abundances of the secondary star. 

In this paper, we use high-resolution spectra to 
derive the stellar parameters and chemical abun- 
dances of the secondary star in the black hole X- 
ray binary V404 Cygni. We then compare in de- 
tail these element abundances in the context of 
the possible enrichment of the secondary star from 
SN/HN yields. 



2. Observations 

We obtained 14 high-quality spectra of V404 Cygni 
in quiescence with the High Re solution Spectrom- 
eter (HIRES; IVoet et al.lll99l at the Keck I 10- 
m telescope (Mauna Kea, Hawaii) on 12 July 
2009 UT. The data covered the spectral regions 
4390-5805 A, 5895-7405 A, and 7520-8780 A 
with a slit width of 0.861", at a resolving power 
A/AA « 50,000. The seeing was in the range 
0.6"-1.0" during the whole night, and for most of 
the main target spectra it was ~ 0.75". This ob- 
serving program was scheduled at a time such that 
the secondary star was near inferior conjunction 
with the black hole (i.e., at an orbital phase of 0), 
to minimize possible effects on the determination 
of chemical abundances due to the asymmetry of 
its Roche lobe. 

We also observed ten template stars with spec- 
tral types in the range K0V/IV-K2V with the 
same instrument and spectral configuration. The 
integration time for V404 Cygni was fixed at 
1800 s in all exposures except for the last two, 
which were of 2100 s and 2000 s duration. This 
relatively long integration time was chosen be- 
cause the orbital smearing at this phase is only 
3-5 kms~^, smaller than the instrumental resolu- 
tion of ^ 6.6 kms~^. 

The spectra were reduced in a standard man- 
ner using the MAKEE package. Each order of 
each individual spectrum was normalized using a 
low-order cubic spline and combined into a single 
one-dimensional spectrum. The individual spec- 
tra were correct ed for their radial velocity (for 



more details, see Gonzalez Hernandez et al. 2011 



in preparation) and combined in order to improve 
the signal-to-noise ratio (S/N). After binning in 
wavelength in steps of 0.1 A, the final spectrum 
had S/N « 70 in the continuum at 6500 A and 
S/N « 150 at 7770 A. This spectrum is displayed 
in Figure [1] in comparison with that of a K2 V 
template star. Note that the spectrum of the sec- 
ondary star in this system and the spectra of the 
template stars were normalized using the same 
procedure. 
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Fig. 1. — Observed spectrum of the secondary star of V404 Cygni, after correcting for the individual radial velocities and 
combining all individual spectra (top panel) , and a properly broadened template also corrected for its radial velocity (HIP 64797, 
bottom panel). 
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Fig. 2. — Distributions obtained for each parameter using Monte Carlo simulations. The bottom-right panel shows the 
distribution obtained for the veiling slope, mo, given as — mo/10~^ in units of A~^. The labels at the top of each bin indicate 
the number of simulations consistent with the bin value. The total number of simulations was 1000. 
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Table 1 

Ranges and Steps of Model Parameters 



Parameter Range Step 



Teff 4200 5800 K 100 K 

log[g/(em s^)] 2^4 0.1 

[Fe/H] -0.5 0.6 0.1 

/4500 0.6 0.05 

mo ^ -0.000171 -0.000019 
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Fig. 3. — Best synthetic spectral fits to the Keck/HIRES spectrum of the secondary star in the V404 Cygni system (bottom 
panel), and the same for a template star (properly broadened) shown for comparison (top panel). Synthetic spectra are computed 
for solar abundances (dashed line) and best-fit abundances (solid line). 



3. Chemical Analysis 
3.1. Stellar Parameters 

The merged and normalized spectrum of the 
secondary star in V404 Cygni may show appar- 
ently weaker stellar lines due to the veiling intro- 
duced by the accretion disk. This veiling is found 
to be ~ 13%±2% in the range 6400-6600 A in a set 
of observations taken duri ng the period 1990-1993 
( Casares fc Charles 19941 ). This veiling was esti- 
mated by performing standard optimal subtrac- 
tion techniques with a KO IV template star. 

The high quality and resolution of the Keck/HIRES 
spectrum enable us to infer the stellar param- 



eters, (Teff, logg, and the metallicity [Fe/H]) 
of the companion star, taking into account any 
possible veiling from the accretion disk as in 
previous studies of other LMXBs (see, e.g., 



IConzalez Hernandez et"al ]l2004l: l iGonzalez Hernandez et al 



.2008bD . For simplicity, the veiling was defined 
as a linear function of wavelength and thus de- 
scribed with two additional parameters, the veil- 
ing at 4500 A, /4500 = J^dfskV^stc™, and the slope. 
Too- We note that the total flux is defined as 
-Ftotai = ^disk + Fscc, where Fdisk and Fscc are the 
flux contributions of the disk and the continuum 
of the secondary star, respectively. 

The most recent version of this code (see 
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Fig. 4. — The same as in Fig. \3\ but for the spectral range 5910-5955 A. In this spectral region there are several relatively 
weak telluric lines labelled as AtmH20. 




Fig. 5. — The same as in Fig.O but for the spectral range 7756-7819 A. 
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also Gonzalez Hernandez et al. I2OO9I) allows us 
to compare, via a ^^"niiiiiniization procedure, 
up to 50 small spectral regions of the stellar 
spectrum with a grid of synthetic spectra com- 
puted using the local thermodynamical equilib- 
rium (LTE) code MOOG (Sncdcn 1973). Wc used 
a gri d of LTE model atmospheres (jKurucz et all 
and the atomic line data from the Vienna 
Atomic Line Database (VALD IPiskunov et al 



19951). The oscillator strengths of relevant lines 
were adjusted until they reproduced the solar at- 



( Grevesse et al 



las (IKurucz et al.l 119841) with solar abundances 
The corrections applied to 



the log 5/ values taken from the VALD database 
were smaller than ~ 0.2 dex. 

We inspected the high-quality HIRES spectrum 
of V404 Cygni, trying to search for Fe I-Fe II lines. 
We finally selected 16 spectral features containing 
more than 50 Fe lines with excitation potentials 
between 1 and 5 eV. We emphasize that these Fe 
lines have a range of log gf values and are spread 
over the whole spectral range 5100-6750 A, and we 
try to find the best-fit model for different veiling 
factors and various veiling slopes. Therefore, the 
strength of a given spectral feature depends not 
only on the stellar parameters and metallicity of 
the star but also on the veiling factor appropriate 
for that feature at the corresponding wavelength. 
In addition, due to the relatively high rotation ve- 
locity of the secondary star, these Fe lines were 
blended with stellar lines of other elements such 
as Ti and Ni. However, the main contributors to 
the selected features were always the Fe lines, and 
the abundances of other elements were in any case 
scaled with the metallicity of the model. 

The five free parameters were varied in the 
ranges given in Table [1] The rotation veloc- 
ity of the companion star was measured to be 
36.4 kms~^ usin g the HIRES spectrum (for fur 



ther details, see Gonzalez Hernandez et al. 2011 



in preparation), and a limb darkening of e = 0.65 
was adopted. The microturbulence, ^, was com- 
puted using an experimental expression as a func- 
tion of effective temperatu re and surface gravity 



(jAUende Prieto et all [20o3 ). The value for the 
best-fit model given below is ^ = 1.117 kms^^. 

We obtain as most likely values T^e = 4800 ± 
100 K, log[g/(cm s2)] = 3.50 ± 0.15, [Fe/H] = 
0.20 ± 0.17, /4500 = 0.15 ± 0.05, and mo = 
-0.00013 ± 0.00002. The la uncertainties of the 



five free parameters were determined using 1000 
realizations whose corresponding histograms are 
displayed in Figure [2] Thus, we find a very small 
veiling, with /5170 ~ 0.05 for the spectral region 
of the Mg I b triplet at 5167-5183 A and almost 
zero for longer wavelengths. This veiling is lower 
than t hat measured during early phases of quies- 



cence (jCasares et al.l 119931 ) . We emphasize that 



the linear function adopted to model the behavior 
of the veiling is strictly valid only in the spectral 
range 5000-6800 A, and therefore should not be 
extrapolated beyond that wavelength range. On 
the other hand, we have analyzed several spectral 
ranges around 7750 A, 8450 A, and 8750 A by as- 
suming a veiling /nir — for all of them, and the 
Fe I lines in these regions seem to be well repro- 
duced by the spectral synthesis. However, a veiling 
of throughout the infrared might not be correct, 
and one would have to analyze spectral lines in 
the infrared to check if this veiling, derived in the 
optical, holds for longer wavelengths. 

The stellar parameters of secondary stars in 
LMXBs are particularly relevant for the determi- 
nation of orbital inclinations, which are typically 
derived from the modelling of ellipsoidal varia- 
tions using opt ical and NIR light curves at quies- 
cence (see, e.g., Gelino et al.ll2006 ). ^Sha hbaz et al. 
(jl994l ) derived Tes ~ 4360 K for the companion 



star in their modelling of the iC-b and light curve of 
V404 Cygni. Hynes e t al.l (|2009[ ) suggest that its 
most likely spectral is KO HI, which implies Toff « 
4570 K using the tem perature sc ale of g i ants i n 
van Belle etall ^M)- iKharghari a et al.l (|2010l ). 
however, propose a spectral type of K3 HI as the 
best template that matches their NIR broad-band 
spectra with little contribution from the accretion 
disk, providing a veiling of ~ 2% and ^ 3% for 
their H— and K-hand spectroscopy. This spec- 
tral type lends to a Tcff ~ 4300 K, according to 
their model. However, our spectroscopic value 
(Tefi « 4800 K) may require more contribution of 
the flux from the accretion disk in the NIR bands 



(see, e.g., Hynes et al. 2005[) . which may also im- 
ply a different orbital inclination, and therefore a 
different black hole mass. 

We have collected the most updated informa- 
tion on orbital parameters for several X-ray bina- 
ries from the literature to derive the expected size 
of the Roche lobe and thus the corresponding sur- 
face gravity. We compare these values with our 
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Table 2 

Chemical Abundances of the Secondary Star in V404 Cygni 



Element 


loge(X)e'' 


[X/H] 


[X/Fc] 


a 


A„ 




'^log g 


As 




A [X/H] 


A[X/Fc] 




Fe 


7.50 


0.23 




0.12 


0.02 


0.05 


0.00 


-0.16 


0.08 


0.19 




51 




8.74 


0.60 


0.37 


0.07 


0.05 


-0.16 


0.06 


-0.06 


0.06 


0.20 


0.24 


2 


Na 


6.33 


0.30 


0.07 


0.20 


0.12 


0.10 


-0.03 


-0.13 


0.10 


0.23 


0.14 


3 


Mg 


7.58 


0.00 


-0.23 


0.09 


0.03 


0.06 


-0.03 


-0.09 


0.08 


0.15 


0.09 


8 


Al 


6.47 


0.38 


0.15 


0.07 


0.05 


0.05 


0.00 


-0.08 


0.05 


0.11 


0.11 


2 


Si 


7.55 


0.36 


0.13 


0.22 


0.06 


-0.06 


0.04 


-0.10 


0.06 


0.15 


0.15 


14 


Ca 


6.36 


0.20 


-0.03 


0.16 


0.05 


0.11 


-0.04 


-0.24 


0.10 


0.29 


0.12 


11 


Ti 


5.02 


0.42 


0.19 


0.24 


0.06 


0.15 


0.02 


-0.28 


0.12 


0.34 


0.17 


15 


Cr 


5.67 


0.31 


0.08 


0.15 


0.07 


0.15 


-0.05 


-0.07 


0.07 


0.20 


0.16 


5 


Ni 


6.25 


0.21 


-0.02 


0.32 


0.11 


0.06 


-0.00 


-0.25 


0.14 


0.31 


0.15 


9 



Note. — Chemical abundances and uncertainties due to the uncertainties At^^.^ — +100 K, Aiogg — 4-0.15 dcx, — 
+0.5 km s~^, and Avci ^ 0.05. 

^T he solar element abun dances were adopted from [Grevesse et al.| jl996D for all elements except oxygen which was taken 
from lEcuvillon et all ( 120061) . 

"^The uncertainties from the dispersion of the best fits to different features, A^r, are estimated using the following formula: 
Act — (t/^/N, where cr is the standard deviation of the measurements. 

'^Number of spectral features of this element analyzed in the star, or if there is only one, its wavelength. 

^The oxygen abun dance has been corrected for NLTE effects, Anlte — —0.2, using the NLTE computations in 
lEcuvillon et al.l ll2006ll . 



Table 3 
Surface Gravities in LMXBs 



Name* 


Porb 

[days] 


Mco.f 
[Mq] 


M2J 

[Mq] 


If 


[RgI 


Rl 

[Rq] 


logs(flL) 


logs(spec) 


V404 Cygni 
Cen X-4 
A0620-00 
Nova Sco 94 
XTE J1118 


6.47129 ± (7 X 10"^) 
0.6290522 ± (4 x 10"'') 
0.32301405 ± (1 X 10"**) 
2.62168 ± (1.4 X 10"*) 
0.1699339 ± (2 X lO"'') 


9.00 ± 0.60 
1.50 ± 0.40 
6.61 ± 0.25 
6.59 ± 0.45 
8.30 ± 0.28 


0.54 ±0.08 
0.30 ± 0.09 
0.40 ± 0.05 
2.76 ± 0.33 
0.22 ± 0.07 


0.060 ± 0.005 
0.200 ± 0.030 
0.060 ± 0.004 
0.419 ± 0.028 
0.027 ± 0.009 


30.98 ± 0.75 
3.76 ± 0.38 
3.79 ± 0.05 

16.85 ± 0.48 
2.64 ± 0.04 


5.51 ± 0.27 
0.95 ± 0.14 
0.67 ± 0.02 
5.17 ± 0.24 
0.37 ± 0.05 


2.69 ± 0.03 
3.96 ± 0.02 
4.38 ± 0.02 
3.45 ± 0.01 
4.65 ± 0.04 


3.50 ± 0.15 
3.90 ± 0.30 
4.20 ± 0.30 
3.70 ± 0.20 
4.60 ± 0.30 



2005), 
2003 ). 



Referenees: V404 Cygn i: Casares ^^Ch^del 994) , Casares (1996), Khargliaria et afl ||20 1^. and ttiis work; Centau r us X- 4: [Torres_e£^ ^ 
_ izalez Hernandez et al.l ||2005) . Casa res et al.|72007), Khar gliaria et al. (2010) : A0620-00: IConzalez Hernandez et al.| p004D.INeilsen et a 
Gonzalez Hernandez fc Casarej pOlO), CantrcU et al. (201 0); Nova S corpii 1994: van dcr Hooft et al. (1998), Gonzalez Hernandez et al. (2008a), 

Torres et al.. 12004) . .Gelino et al.l I I2006I) . [Gonzalez Hernandez et al.. (.2006.) . .Gonzalez Hernandez et al.. (.2008h) 



Shahbaz (2003); XTE J1118+480: 
and this work; 



Note. — Comparison between LMXB surface gravity determined spectroscopically, log (/(spec), with the surface gravity, log g^R.^) , computed from 
the orbital period, Porb, the mass ratio, q/, and the current masses of the compact object, iV/co,/, and the secondary star, A^2,/. 
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Table 4 

kinematical properties in lmxbs 



Name* 


K2 




/(M) 


V sin i 


7 


d 




Mi 






[km 




[Mq] 


[kms^i] 


[kms^i] 


[Kpc] 


[mas yr~ ^] 


[mas yr 


-'\ 


V404 Cygni 


208.4 ± 


0.6 


6.09 ± 0.04 


40.8 ± 0.9 


0.3 ± 0.6 


2.39 ± 0.14 


-5.04 ± 0.02 


-7.64 ± 


0.03 


Cen X-4 


144.6 ± 


0.3 


0.197 ± 0.001 


44 ± 3 


189.6 ± 0.2 


1.30 ± 0.40 


11 ± 10 


-56 ± 


10 


A0620-00 


435.4 ± 


0.5 


3.10 ± 0.04 


82 ± 2 


8.5 ± 1.8 


1.06 ± 0.12 








Nova Sco 94 


226.1 ± 


0.8 


3.14 ± 0.03 


86 ± 4 


-141.9 ± 1.3 


3.20 ± 0.50 


-3.3 ± 0.5 


-4.0 ± 


0.4 


XTE J1118 


708.8 ± 


1.4 


6.27 ± 0.04 




2.7± 1.1 


1.70 ± 0.10 


-16.8 ± 1.6 


-7.4 ± 


1.6 



* Referenc es: V 404 Cygni: 'CasarcS ('1996'), 'Casarei (12007^, 'Miller- Jones et al.l |I2009H), 'Khargharia et al.' f2010'); Cen- 
taurus X-4: IGon zalez Hernandez et al. (2005), Casarcs ct al. (2007), Kliargharia et al. (2010); A0620-00: Neilsen et 
1200811, 'Gonzalez Hernandez Casarcs (2010), Cantrell et al. (2010); Nov a Seor pii 1994: iHiellming &: Rupen (,1995ir 
Orosz fc Bailvn (1997), Mirabel et al. (2002), Gonzalez Hernandez et al. I2008al ). and this work; XTE J1118-I-480: 
Mirabel et al.. (.2001' ), .Gelino et al., (,2006i) . iGonzalez Hernandez et al.i t2008hi) . 



Note. — Kinematical and dynamical properties, orbital parameters and proper motions of low-mass X-ray binaries. 



Table 5 

Stellar, veiling parameters and chemical abundances in LMXBs 



Name* A0620-00 Centaurus X-4 XTE ,71118+480 Nova Sco 94 V404 Cygni 



Alternative name 


V616 Mon 


V822 Cen 




GRO J1655-40 


GS 2023-^338 


Teff (K) 


4900 ± 100 


4500 ± 100 


4700 ± 100 


6100 ± 200 


4800 ± 100 


log(9/cm s^) 


4,2 ± 0,3 


3,9 ± 0,3 


4,6 ± 0,3 


3,7 ± 0,2 


3,50 ± 0,15 


/4500 


0,25 ± 0,05 


1,85 ± 0,10 


0,85 ± 0,20 


0,15 ± 0,05 


0,15 ± 0,05 


mo/10~* 


-1,4 ± 0,2 


-7,1 ± 0,3 


-2 ± 1 


-1,2 ± 0,3 


-1,3 ± 0,2 


[0/H]t 








0,91 ± 0,09 


0,60 ± 0,19 


[Na/H] 








0,31 ± 0,26 


0,30 ± 0,19 


[Mg/H] 


0,40 ± 0,16 


0,35 ± 0,17 


0,35 ± 0,25 


0,48 ± 0,15 


0,00 ± 0,11 


[Al/H] 


0,40 ± 0,12 


0,30 ± 0,17 


0,60 ± 0,20 


0,05 ± 0,18 


0,38 ± 0,09 


[Si/H] 






0,37 ± 0,21 


0,58 ± 0,08 


0,36 ± 0,11 


[S/H] 








0,66 ± 0,12 




[Ca/H] 


0,10 ± 0,20 


0,21 ± 0,17 


0,15 ± 0,23 


-0,02 ± 0,14 


0,20 ± 0,16 


[Ti/H] 


0,37 ± 0,23 


0,40 ± 0,17 


0,32 ± 0,26 


0,27 ± 0,22 


0,42 ± 0,20 


[Cr/H] 










0,31 ± 0,19 


[Fe/H] 


0,14 ± 0,20 


0,23 ± 0,10 


0,18 ± 0,17 


-0,11 ± 0,10 


0,23 ± 0,09 


[Ni/H] 


0,27 ± 0,10 


0,35 ± 0,10 


0,30 ± 0,21 


0,00 ± 0,21 


0,21 ± 0,19 



References: V404 Cygni: This work; Centaurus X-4: 



pOOTI).: A0620-00: |Gonzalez Hernandez et all (120041), 


Gonzalez Hernandez et al. 


IConzalez Hernandez et all ll2008a|l: XTE J1118+480: 


Conzalez Hernandez et al. 



^Thc uncertainties on the stellar abundances given in this table has been derived without taking into the 
error on the microturbulence (sec Table [21 . 

^Oxygen abundances are given in NLTE. 



Note. — Stellar and veiling parameters, and chemical abundances of secondary stars in low-mass X-ray 
binaries. 
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spectroscopic determinations in Table [5l Using 
Kepler's third law, we have estimated the current 
orbital separation, Oc,/, from the orbital period, 
Porb, and the reported masses of the compact ob- 
ject, Mcoj, and the companion star, M2J. The 
ratio of the radius of the Roche lobe, Rl, and the 
orbital separatio n can be estini ated using Eggle- 
ton's expression ( Eggletonl 1983 ) as 



RL/acj = 0.49g^/7[0.6g^/^ + In (1 + q''% 

Assuming that the star is filling its Roche lobe, 
Rl = R2J, and using the mass of the secondary 
star, we derive the expected surface gravity. 

From Table [5] one can see that in general, the 
spectroscopic determinations are close to values 
derived from the size of the Roche lobe of the 
secondary star. In particular, for Centaurus X-4, 
A0620-00, and XTE J1118+480, the values are 
consistent within the la uncertainties, and for 
Nova Sco 94 at 1.2a. However, for V404 Cygni, 
the surface gravity is only compatible at the 5a 
level, given the relatively small uncertainty in the 
spectroscopic surface gravity. On the other hand, 
we note that surface gravities derived from the 
size of the Roche lobe strongly depend on the cur- 
rent masses of the compact object and the sec- 
ondary star, and thus depend on the estimated, 
sometimes uncertain, orbital inclination. We note, 
for instance, the case of A0620-00, whose black 
hole mass estimate has ch anged over the last ten 
years from 11.0 ± 1.9 M,7. (ICelino et al.ll200ll) . to 



9.7 ± 0.6 Mq (iFroning et al.l 120071). an d finally 
to 6.61 ± 0.25 MrT^ (|Cantren et al.l l2010l) . How- 
ever, for the case of V404 Cygni, the mass ra- 
tio and the orbital period have been determined 
with high precision (see Table [5]), leaving little 
room to increase the value of the surface grav- 
ity obtained from the secondary's mass and the 

size of the Roche lob e. Thus, adopting larger 

masses ( Shahbaz et al. I ll994h for the compact ob- 
ject, Mcoj = 12 Mq, and the companion star, 
M2J = 0.7 Mq, the result does not change sig- 
nificantly; the estimated surface gravity would be 
loggRi = 2.72 ± 0.13. Therefore, a possible un- 
certain orbital inclination does not seem to be the 
reason for this disagreement. 

3.2. Stellar Abundances 

We inspected several regions in the observed 
Keck/HIRES spectrum of the secondary star, 



searching for suitable lines for a detailed chem- 
ical analysis. Using the derived stellar param- 
eters, we first determined the Fe abundance by 
comparing synthetic spectra with each individ- 
ual feature in the HIRES spectrum (see Table [2]). 
In Figure [3] we display one of the spectral re- 
gions analyzed to obtain the Fe abundance, also 
showing the best synthetic spectral fit to the ob- 
served spectrum of a template star (HIP 64797 
with Toff^ 4970 K, log g= 4.33, and [Fe/H] 
= —0.24 dex). We only use as abundance indi- 
cators those features which are well reproduced 
in the template star. The chemical analysis is 
summarized in Table [2] The errors in the element 
abundances show their sensitivity to the uncer- 
tainties in the effective temperature ( At^jj. ) , grav- 
ity (Aiogg), veiling (Avoi), microturbulence (A^), 
and the dispersion of the measurements from dif- 
ferent spectral features (Ao-). In Table [2] we also 
state the number of features analyzed for each el- 
ement. The uncertainties were estimated as 
Act — a/^/N, where a is the standard deviation 
of the N measurements. The uncertainties Aj-^fj , 
A^, Aiogg, and Avoi were determined in the same 
way as, for instance, in the Tcff case: Aj-^jj = 

{J2i=i ATett,i)/-^- The total uncertainty given in 
Table[2]was derived using the following expression: 

A[X/H] = ^^7^~T^~^^7^,. 

In Figure m we show the spectral region 5910- 
5955 A, where there are some Ti and Si lines used 
to derive their abundances. In this region there 
are also some telluric lines, with equivalent widths 
of ~ 30 mA, which are significantly weaker than 
the stellar features, whose equivalent widths are 
^ 250 mA, and we do not think these lines have 
caused any problem in the abundance determina- 
tion. This spectral region contains two Fe features 
at ^ 5914.2 and '-^ 5916.3 A which are sensitive 
to surface gravity variations. One can easily no- 
tice the relative different strengths of both Fe fea- 
tures in the template spectra (with log g« 4.3) 
and in the secondary star in V404 Cygni (with 
log gfn 3.5). 

3.3. Oxygen 

The oxygen abundance was derived from the 
O I triplet at 7771-5 A, which for the relatively 
high rotation velocity of the secondary star in 
V404 Cygni produces two well-resolved features 
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that we have analyzed independently. In Figure [5] 
we show the 7755-7820 A range, where the O I 
features in the secondary star appear to be appar- 
ently enhanced when compared with those in the 
template star having similar effective temperature. 
The best fit in LTE gives an oxygen abundance of 
[0/H]lte ~ 0.8. We note that the atomic data 



for O I lines were adopted from lEcuvillon et al 
( 2006h . These authors slightly modified the oscil- 



lator strengths in order to obtain a solar oxygen 
abundance of loge(O)0 = 8.74. For other ele- 
m ents we assumed as sola r abundances those given 
bv lGrevesse et al.l (|l996[ ). 

The O I A7771-5 triplet suffers from ap- 
preciable non-LTE (NLTE) effects (see, e.g., 



Ecuvillon et al. 20061) . For the stellar param- 
eters and oxygen abundance of the secondary 
star, NLTE correctiontQ are estimated to be 

0.20 dex for the NIR O I A7771-5 triplet. 

Table [5] provides the oxygen abundance properly 
corrected for NLTE effects. 

In principle, one could argue that the relatively 
"high" value for the spectroscopic surface grav- 
ity found in this work (see i 33.ip for V404 Cygni 
may affect the derived abundances. We note that 
the spectroscopic values were derived for the sec- 
ondary star at nearly inferior conjunction, which 
means that we are looking at the "spherical" side 
of the Roche-lobe-like star. Therefore, the derived 
Tcs and log g values should be larger than the 
mean Toff and log g values of the star (see Ta- 
ble[5]), although this would only account for 100 K 
and 0.1 dex, respectively. 

In Table [2] one can see that the NIR O I lines 
are sensitive to the surface gravity, with a change 
of 0.06 dex in derived abundance for a change of 
-1-0.15 dex in \ogg. This means that for a surface 
gravity as low as \ogg = 2.7, one would obtain a 
0.32 dex lower oxygen abundance. On the other 
hand, the Fe abundance would not be sensitive to 
this change in surface gravity, and therefore we 
would get [0/Fe]NLTE = 0.05 dexH We also note 
that, for instance, a 300 K lower Tgg would re- 



^AnlTE = loge(X)NLTE - log<;(X)LTE- 

■^Fe lines are indeed quite sensitive to variations of the sur- 
face gravity, typically with changes of ±0.05 dex in Fe abun- 
dance for a change of -1-0.15 dex in log 3, although some Fe 
lines would not show any change in abundance. However, 
these differences compensate each other, so that the aver- 
age change is ~ 0.001 dex, almost null. 



suit in a substantial increase in the derived oxygen 
abundance, ~ 0.48 dex (see Table [2]), and a smaller 
decrease in the Fe abundance, ~ —0.15 dex, yield- 
ing [0/Fe]NLTE = 1 dex with [Fe/H] = 0.08. 
These values has been estimated without taking 
into account the possibly smaller NLTE correction 
at the lower value of Toff. 

3.4. Magnesium 

The magnesium abundance was derived from 
eight Mg I features, including the four opti- 
cal lines of the Mg I b A5167-83 A triplet and 
Mg I A5528 A. These Mg I lines ar e known to 
be se nsitive to NLTE effects (see, e.g., Zhao et all 
19981 ). although with NLTE corrections of only 



~ -1-0.05 in the Sun and p robably smaller ones a t 
these cooler temperatures ( Zhao fc GehrenI 2000[ ). 

We note here that we are applying an auto- 
matic 1.5(7 abundance rejection over the initial set 
of lines, slightly affecting the Mg abundance. This 
discards two Mg I features, lowering the Mg abun- 
dance by 0.06 dex and decreasing the standard de- 
viation of the measurements from 0.15 to 0.09 dex. 
For other elements like Si and Ti, the 1.5cr rejec- 
tion yields 0.04 lower and 0.03 higher abundance, 
respectively, with a decrease of the standard devi- 
ation by only 0.02 dex. 

Finally, possible differences in the stellar pa- 
rameters also have an impact on the Mg abun- 
dance (see Table[2]). The Mg I lines are sensitive to 
the effective temperature and surface gravity, with 
a change of 0.06 and —0.03 dex in derived abun- 
dance for changes of 100 K in Toff and -1-0.15 dex 
in log g, respectively. In this case, 300 K lower 
Toff would result in 0.18 dex lower Mg abundance, 
whereas a decrease of 0.8 dex in log g would give 
rise to a 0.16 dex higher Mg abundance. 

4. Discussion 

The global metallicity of the secondary star in 
V404 Cygni is slightly higher than solar, similar to 
that of the secondary star in other black hole X-ray 
binaries such as A0620-00 and XTE Jl 118-1-480, 
but also comparable to that of many stars in the 
solar neighborhood. We have searched for anoma- 
lies in the abundance pattern of the secondary star 
in comparison with typical abundances of stars be- 
longing to the Galactic thin disk. The abundances 
of other elements relative to iron listed in Table [2] 
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Fig. 6. — Abundance ratios of tlie secondary star in V404 Cygni (blue wide cross) in comparison with the abundances of 
solar-type metal-rich dwarf stars. Galactic trends were taken from [Gonzalez Hernandez et all l l2010t) . The size of the cross 
indicates the uncertainty. Filled and empty circles correspond to abundances for exoplanet host stars and stars without known 
exoplanet companions, respectively. The dashed-dotted lines indicate solar abundance values. 



are compared in Figures |6] and [7] with the Galactic 
trends of these elements in the relevant range of 
metallicities. We ado pted recent, very accurate 



Galactic trends from Gonzalez Hernandez et al 



( 20101 ) for all elements except for oxygen, because 
we needed O abundances measured from the NIR 
O I triplet. We thus decided to a dopt the Galac- 
tic tre nd of oxygen published by lEcuvillon et al 



pood ). Figures [B] and [7] show moderate anoma- 
lies for Ti, Al, and Si in the secondary star in 
V404 Cygni, and whereas Mg is surprisingly un- 
der abundant, oxygen appears to be enhanced in 
comparison with the galactic trend. 

Table H] shows the element abundance ratios in 
V404 Cygni and the average values in stars with 
iron content in the range 0.11 < [Fe/H] < 0.35, the 
comparison sample corresponding to a la uncer- 
tainty in the iron abundance of the secondary star. 
The elements Na, Si, Ca, Cr, and Ni are consistent 
with the average values in thin-disk stars within 
the error bars. Mg appears to be under abundant, 
whereas Al and Ti are roughly consistent at Icr. 



Oxygen may be considered separately since it is 
clearly more abundant, at 2.8cr, than the average 
values of the thin-disk stars. 
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Fig. 7. — The same a s in Fig. [6] but for oxyg en. Galactic 
trends were taken from lEcuvillon et al.l ||2006| ). 



V404 Cygni is located in a Galactic thin-disk 
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Table 6 



Element Abundance Ratios in V404 Cygni 


Element 


[X/Fe]„o,cyg 


^[X/Fo],V404Cyg 








O 


0.37 


0.24 


-0.32 


0.19 


0.033 


Na 


0.07 


0.14 


0.05 


0.05 


0.009 


Mg 


-0.23 


0.09 


0.00 


0.05 


0.008 


Al 


0.15 


0.11 


0.04 


0.03 


0.006 


Si 


0.13 


0.15 


0.00 


0.01 


0.002 


Ca 


-0.03 


0.12 


-0.02 


0.03 


0.005 


Ti 


0.19 


0.17 


0.01 


0.02 


0.003 


Cr 


0.08 


0.16 


0.00 


0.01 


0.001 


Ni 


-0.02 


0.15 


0.03 


0.02 


0.003 



* Uncertainties in the element abundance ratios ([X/Fe]) in the secondary star in 
V404 Cygni. 



Note. — [^/Fe]^^^^^^ indicates the average value of 36 stars with iron content in the 
range 0.11 to 0.35 corresponding to the abundance range [Fe/H]±(TFc of the secondary 
star in V404 Cygni, taken from Ecuvillon et al. ( 2006) and Gonzalez Hernandez et al.| 
l|2010() . The uncertainty in the average value of an abundance ratio in the comparison 
sample is obtained as A^-, stars — "^stars/V^, where crgtars is the standard deviation 
of the measurements and A'^ is the number of stars. 



regio n at a distance d = 2.39 =b 0.14 kpc from the 
Sun (jMiller- Jones et al.l |2009b|) , but its peculiar 



velocity is larger than the int rinsic velocity dis- 
persion in the Galactic plane. iMiller -Jones et al 



argue that this pecuhar velocity, I'pec — 
39.9 ± 5.5 kms~^, may be the result of the na- 
tal kick from the mass ejection in the SN/HN ex- 
plosion of the primary star in this LMXB. The 
companion star could have captured a significant 
amount of the ejected matter in the SN/HN ex- 
plosion that formed the compact object. We have 
explored this possibility with the aim of obtaining 
information on the chemical composition of the 
progenitor of the compact object. 

4.1. Spherical Explosion 

The secondary star in V404 Cygni may have 
lost a significant amount of its initial mass 
through a mass-transfer mechanism onto the 
compact object, during the binary evolution. 
Its post -SN evolution has been studied in de- 
tail by iMiller- Jones et al.l (|2009al) . These au- 
thors suggest that from 0.5 to 1.5 M© may have 
been accreted from the companion star. We 
estimated the current mass of the secondary 
star at M2J = 0.54 ± 0.08 M© (see Table EJ, 
from the cur rent black hole mass Mbh,/ = 
9.0 ± 0.6 Mq (jKhargharia et al.l I2OIOI) and mass 
ratio qf = 0.060 zb 0.005 (jCasares fc Charles! 



19941 ). Following IMiller- Jones et al.l (l2009a|), we 
assume an initial mass for the secondary star of 
M2,i « 2 M0 and an initial black hole mass of 
A^BH.i ~ 8 Mq. Therefore, the final mass of the 
black hole would be Mbhj ~ 9.5 Mq which is 
consistent with the current estima te of the black 
hole mass ( Khargharia et al. 2010l ). 



The high O content in this companion star may 
not suggest any strong CNO processing within the 
star itself during its evolution, since this w ould in- 
creas e the N abundance (see, e.g., .Has well et al 
20021 ). Oxygen, however, is not expected to de - 



crease too much compared withC ( Clayton 19831) . 
Therefore, the CNO-processed material is C un- 
derabundant and N overabundant. Unfortunately, 
in the Keck/HIRES optical spectrum of the sec- 
ondary star there are no sufficiently strong and 
unblended stellar lines t o provide an accur a te and 



reliable C abundance. iKhargharia et al. 



argued that their NIR spectrum of V404 Cygni 
suggests that CO molecules seems to fit well with 
a solar-abundance K3 III template spectrum for 
which they adopt Tcs « 4300 K. We note here 
that CO molecular bands are more sensitive to 
the C abundance than to the O abundance and 
are stronger at lower effective temperatures. Thus, 
their argument do not necessarily contradict our 
determination of the O abundance. 

The present orbital distance is Ucj ~ 31 R© (see 
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Table [5]). We will assume that the post-SN orbital 
separation after tidal circularization of the or- 
bit was ac,i ~ 15 R©, since the secondary must 
have experienced mass and angular momentum 
losses during the binary evolution until it reached 
its present configurati on. Similarly to the case 
of XT E Jl 118+480 (|Gonzalez Hernandez et"aD 
2008bl ). we can estimate the maximum ejected 
mass in the SN/HN explosion. A binary system 
such as V404 Cygni will survive a spherical SN ex- 
plosion if the eject ed mass AM = Mne — Mbh,! < 
(A/hc + M2,0/2 (|Hillsl Il983[ ). This implies a 



mass of the He core before the SN explosion of 
< 18 Mq. We therefore adopt He core 
masses of Mho ~ 15-16 M©, so we infer a He 
co re radius of i?He ~ 3 R© from the expression 



Portegies Zwart et al.l (|1997al . and references 



therein). 

Assuming a pre-SN circular orbit and an instan- 
taneous spherically symmetric ejection (that is, 
shorter than the orbital period), one can estimate 
the pre-SN orb ital separation, an, using the re - 
lation given by Ivan den Heuvel fc HabetsI (Il984h : 
ao = ac,i lJ.f, where fj,f = (Mbh,! + -M'2,i)/(MHe + 
Af2,i)- We find flQ « 9, for the adopted values 
of Mhc = 15 Mq and a^i = 1 5 R©. At the time 
of the SN explosion (~5-6 Mvr: lBrunish fc Truran 
19821 ). a 2 Mq secondary star, still in its pre-main- 



sequence evolution, has a radius i?2,i ~ 3.3 R©. 
The material possibly captured by the secondary 
star in the SN/HN event has a much larger mean 
molecular weight than the pre-explosion gas of the 
secondary star, so it has probably been well mixed 
within the star by ther mohaline mixing in a rela- 



tively short timescale (jPodsiadlowski et al.l 12002 . 



and references therein). Therefore, the subse- 
quent lost of material from the secondary star onto 
the compact object should not affect the chemi- 
cal composition of its atmosphere. The amount of 
mass deposited on the secondary can be estimated 

as ^Tlcap 

= AM{TrRlJ4Tral)fce.p M©, where /cap 
is the fraction of mass, ejected within the solid an- 
gle subtended by the secondary star, that is even- 
tually captured. We assume that the captured 
mass, rricap, is completely mixed with the rest of 
the companion star. 

We compute the expected abundances in 
the atmosphere of the secondary star after the 
pollution fro m the progenitor of the c o mpact 
object as in Gonzalez Hernandez et al. ( 20041 ) 



and Gonzalez Hernandez et al. I (HooiE). We use 
40 M© spherically symmetric core-collapse explo- 
sion models (Mho « 15.1 M©) at solar metallicity 
(Z = 0.02) and for two dif f erent explosion energies 
( Umeda fc Nomotol l2002i 120051: iTominaga etall 
2007t ). These models imply AM 7 M© and 



need small capture efficiencies of /cap ^0.1 (i.e., 
10%) to increase significantly the metal content 
of the secondary star. In order to fix the /cap 
parameter, we have tried to get an expected Al 
abundance, [Al/H] « 0.45, consistent with the ob- 
served Al abundance within the error bars. These 
model computations are shown in Table [T] These 
models would also provide a different mass frac- 
tion of each element at each value of the mass 
cut (defined as the mass that initially collapsed to 
form the comp act remnant). For mor e details of 
the models, see Tominag a et all (|2007h . 

The assumptions regarding the initial mass and 
the post-SN orbital distance of the secondary star 
are not so relevant due to the free parameter /cap- 
A different value for the initial mass of the sec- 
ondary star, as small as (say) 1 M© with an initial 
radius of i?2.i ~ 1-3 R©at the time of the SN ex- 
plosion, and for the post-SN orbital separation, 
say ac,i = 25 R©, would require a larger capture 
efficiency factor of /cap « 0.45 (i.e., 45%). 

The explosion energy is Ek = 1 x 10^^ 
erg and Ek = 30 x 10^^ erg for the SN and 
HN models, respectively. This energy is de- 
posited instantaneously in the central region of 
the progenitor core to generate a strong shock 
wave. The subsequent propagation of the shock 
wave is followed throug: h a hydrodynamic code 
(lUmeda. fc N omoto 2 002[ and refer e nces th erein). 
As in Gonzalez Hernandez et al. I (l2008bh . our 



model computations assume different mass cuts 
and fallback masses, and a mixing factor of 1 
which assumes that all fallback matter is well 
mixed with the ejecta. The amount of fallback, 
Mfaii, is the difference between the final remnant 
mass, Mbh,!, and the initial remnant mass of the 
explosion, Mcut- We should recall here the ejected 
mass, AM, which is equal to Mhc — Mbh.i, where 
Mhc is the mass of the He core. 

We use SN/HN models to provide us with the 
yields of the explosion before radiative decay of 
element species. We then compute the integrated, 
decayed yields of the ejecta by adopting a mass 
cut and by mixing all of the material above the 
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mass cut. Finally, we calculate the composition 
of the matter captured by the secondary star, 
and we mix it with the material of its convec- 
tive envelope. It has been suggested that the 
black hole in the LMXB Nova Sco 1994 could have 
formed in a two-stage process where the initial 
collapse led to the formation of a neutron star 
accompanied by a substantial kick and the fi- 
nal mass of the compact remnant was achieved 
by m atter that fell back after the initial col- 
lapse (jPodsiadlowski et al.l |2002|) . However, the 



black hole mass in that system has been esti- 
mated at ~ 5. 4 Mr7^ dfieer fc Podsiadlowskill2002h 
or ~ 6.6 Mq (|Shahbajl2003l . see Table [5|). while 
the black hole in V404 Cygni may have an initial 
mass of ^ 8 Mq which would require a fallback 
mass of ~ 6.6 Mq if we assume ~ 1.4 M f^ for a 



canonical neutron star. iMacFadven et al.l (|200ll ) 

proposed a scenario where collapsar models har- 
bor black holes that could form in a mild explosion 
with substantial fallback (up to '-^ 5 Mq ). 

In Figure |8] we show the expected abundances 
of the secondary star after contamination from 
the nucleosynthetic products of the SN explosion 
{Ek — Ix 10^^ erg) of a A/hc ~ 15 Mq progenitor 
star. The initial abundances of the secondary star 
have been estimated from the average abundances 
of thin-disk solar-type stars wit h [Fe/H] 0.23 
( Gonzalez Hernandez et al. 2010l ). which are pro- 
vided in Table [T] Note that for a given model, 
the Al abundance in the secondary star hardly de- 
pends on the mass cut, since Al forms in the outer 
layers of the explosion. Thus, once the capture 
efficiency is fixed, the Al abundance in Figure [8] is 
almost constant. 

The expected abundances of O, Mg, and Al 
in the secondary star after contamination from 
the SN ejecta remain mostly independent of the 
adopted mass cut, whereas other elements like Si, 
Ca, Ti, Fe, and Ni are quite sensitive to the mass 
cut of the model. This appears to be more clear in 
the HN model, where the higher explosion energy 
enhances the amounts of the a-process elements 
Si, Ca, and Ti. The Fe-group elements like Cr, 
Fe, and Ni are also sensitive to the mass cut, es- 
pecially in the HN model. 

For both the SN model (left panel in Fig. H]) 
and the HN model (right panel in Fig.[8l), the pre- 
dicted abundances agree with the observed abun- 
dances in the secondary star relatively well, except 



for those of O and Mg. The expected Mg abun- 
dance is too high in comparison with the observed 
value since the adopted initial Mg abundance is 
already higher than the observed abundance (see 
Table [7]). The case of O is just the contrary; the 
initial abundance is so low that it is not possible 
to reach the observed O abundance when fitting 
other element abundances (see Fig. |8|). 

For mass cuts above ~ 3 Mq, very little Si, Ca, 
Ti, Fe, Cr, and Ni is ejected; thus, the expected 
abundances of the model essentially reflect the ini- 
tial abundances of the secondary star. In contrast, 
O, Na, Mg, and Al are hardly sensitive to the mass 
cut, and all are enhanced due to the capture of en- 
riched material in the SN/HN explosion. 

We could have also investigated a model with 
solar initial abundances for the secondary star 
([X/H]g = 0) and the same explosion model of so- 
lar metallicity {Z = 0.02), although it is unlikely 
that a significant amount of elements formed in 
the inner layers of the explosion (such as Ti, Ni, 
and Fe) would be present in the ejecta. In such 
models, only the HN explosion with a mass cut 
as low as Mcut ~ 2 M© would fit the relatively 
enhanced observed abundances of Si, Ti, Cr, Fe, 
and Ni, but we would still not be able to fit the 
O and Mg abundances. This low mass cut would 
require very efficient mixing processes, and the 
model would also need a small capture efficiency 
of /cap ~ 0.15 (i.e., only 15% of the matter ejected 
within the solid angle subtended by the secondary 
star). 

4.2. Aspherical Explosion 

Miller- Jones et al. pointed out that due 



to the relatively small peculiar velocity of the sys- 
tem, V404 Cygni does not require an asymmet- 
ric kick. However, here we explor e this possi- 



bility using explosion models from iMaeda et al 



(j2002[ ) that are not spherically symmetric. An 
aspherical SN explosion produces chemical inho- 
mogeneities in the ejecta which are dependent on 
direction. Thus, if the jet in the aspherical SN 
explosion is collimated perpendicular to the or- 
bital plane of the binary (for more details, see 
IConzalez Hernandez et aL 2005), where the sec- 
ondary star is located, elements such as Ti, Fe, 
and Ni are ejected mainly in the jet direction, 
while O, Mg, Al, Si, and S are preferentially 
ejected near the equatorial plane of the helium star 
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Table 7 

Metal- Rich Supernova/Hypernova Explosion Models in V404 Cygni 



[X/H] Expected'' 

Element [X/H] Observed'' [X/H](,'' Mem.iow" M<,„t,up M,„t,io„ M<,„t, 

Spherical explosion model of Z = 0.02 
Supernova Hypernova 



o 


0.60 


-0.09 


0.13 


0.12 


0.13 


0.12 


Na 


0.30 


0.29 


0.56 


0.56 


0.50 


0.56 


Mg 


0.00 


0.24 


0.35 


0.34 


0.35 


0.34 


Al 


0.38 


0.27 


0.46 


0.45 


0.44 


0.45 


Si 


0.36 


0.24 


0.30 


0.25 


0.34 


0.25 


Ca 


0.20 


0.22 


0.26 


0.22 


0.29 


0.22 


Ti 


0.42 


0.25 


0.28 


0.25 


0.35 


0.25 


Cr 


0.31 


0.24 


0.28 


0.24 


0.31 


0.24 


Fe 


0.23 


0.23 


0.32 


0.23 


0.34 


0.23 


Ni 


0.21 


0.27 


0.41 


0.28 


0.39 


0.28 



Aspherical explosion model with Z = 0.02 
Angle"* = 0°-15° Angle = 0°-90° 



o 


0.60 


-0.09 


0.17 


0.14 


0.19 


0.19 


Na 


0.30 


0.29 


0.39 


0.38 


0.38 


0.39 


Mg 


0.00 


0.24 


0.36 


0.35 


0.36 


0.37 


Al 


0.38 


0.27 


0.45 


0.43 


0.45 


0.46 


Si 


0.36 


0.24 


0.31 


0.28 


0.39 


0.33 


Ca 


0.20 


0.22 


0.24 


0.22 


0.33 


0.27 


Ti 


0.42 


0.25 


0.25 


0.25 


0.41 


0.35 


Cr 


0.31 


0.24 


0.24 


0.24 


0.33 


0.29 


Fe 


0.23 


0.23 


0.23 


0.23 


0.31 


0.26 


Ni 


0.21 


0.27 


0.27 


0.27 


0.40 


0.33 



^Observed abundances of the secondary star in V404 Cygni. 

^Initial abundances assumed for the secondary star in V404 Cygni; see text. 

■^Mcut.iow and Mcut.up are (respectively) the lower and upper mass cuts adopted in 
the model computation. See the exact value in the captions of Figs. 12 and 13. 

■^Expected abundances of the secondary star. 

^Angular range, measured from the equatorial plane, in which all the ejected material 
in the explosion has been completely mixed for each velocity point. 

Note. — Expected abundances in the secondary atmosphere contaminated with nu- 
cleosynthetic products of metal-rich explosion models for two different mass cuts and 
symmetries, presented in Figs. 12 and 13. 
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Fig. 8. — Left panel: Observed abundances (filled circles witii error bars) in comparison with the expected abundances in the 
secondary star of V404 Cygni after having captured 10% of the matter ejected within the solid angle subtended by the secondary 
from a solar metallicity {Z = 0.02) 40 Mq spherically symmetric supernova explosion (Mjje = 15.1 Mq) with Ek = 1 x 10^^ erg 
for two different mass cuts, Mcut = 1.46 Mq (solid line with open circles) and Mcut = 8.09 Mq (dashed-dotted line with open 
circles) . The initial abundances of the secondary star were adopted for the average abundances of thin-disk solar- type stars with 
[Fe/H] = 0.23 ± 0.12. Right panel: Same as left panel, but for a spherically symmetric hypernova explosion (Ek = 30 x 10^^ 
erg) for two different mass cuts, Mcut = 1.74 Mq (solid line with open circles) and Mcut = 8.09 Mq (dashed-dotted line with 
open circles). 



(jMaeda et al.ll2002l ). 

In Figure [S] we compare the predicted abun- 
dances in the atmosphere of the secondary star af- 
ter pollution from an aspherical explosion model 
of a metal-rich progenitor having a 16 Mq He 
core (see also Table [7|). The initial abundances 
of the secondary were, as in the spherical case, 
equally extracted from the average values of solar- 
type stars of the solar neighborhood with simi- 
lar iron content (see Table [7]). The left panel of 
Figure |9] reflects the composition of the material 
ejected in the equatorial plane, which should be 
rich in O, Na, Mg, Al, and Si — the only elements 
significantly enhanced with respect to the initial 
abundances. We see that the choice of the mass 
cut does not make a significant difference, and as 
in the spherical case, this model fits all of the el- 
ement abundances within their error bars except 
for O and Mg. 

In the right panel of Figure |9] we have c onsid- 
ered complete lateral mixing (Podsiadlow ski et al.l 
20021 ) — that is, the ejected matter is completely 
mix ed within each velocity bin (for m ore details, 



see 



Gonzalez Hernandez et all 2008b[) . The ob 



served abundances might be better reproduced if 
complete lateral mixing is adopted, since this pro- 
cess tends to enhance the Si, Ca, Ti, Cr, Fe, and 



Ni element abundances at all mass cuts. A model 
with solar initial abundances for the secondary 
star (i.e., [X/H]p = 0) and the same explosion 
model of solar metallicity {Z = 0.02) was also in- 
spected, providing the same result except for the 
equatorial model, which produces too low abun- 
dances of Ti, Fe, and Ni in comparison with the 
observations. 

5. Conclusions 

We have presented Keck I/HIRES high- resolution 
spectroscopy of the black hole X-ray binary 
V404 Cygni. The spectra were obtained at system 
orbital phase close to zero to minimize the effect 
of the Roche-lobe symmetry of the secondary star. 
We have performed a detailed chemical analysis 
of the secondary star, applying a technique that 
provides a determination of the stellar parameters 
and takes into account any possible veiling from 
the accretion disk. We find Tes = 4800 ± 100 K, 
log[g/cm s2] = 3.50 ± 0.15, [Fe/H] = 0.23 ± 0.19, 
and a disk veihng (defined as Fdisk/^totai) of ^ 5% 
at 5100 A, decreasing toward longer wavelengths. 

We have derived the chemical abundances of 
O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, and Ni. They 
are typically higher than solar, and some elements 
show additional slight enhancements (e.g., Al, Si, 
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Fig. 9. — Left panel: Observed abundances (filled circles with error bars) in comparison with the expected abundances in 
the secondary star of V404 Cygni after having captured ~ 10% of the matter ejected within the solid angle subtended by 
the secondary from an aspherical SN explosion of Ek = 10 x 10^^ erg for two different mass cuts, Mcut = 2.41 Mq (solid 
line with open circles) and Mcut = 7.56 Mq (dashed-dotted line with open circles). This model corresponds to the matter 
ejected in the equatorial plane of the primary, where we assumed that the secondary star is located (for more details, see 
iGonzalez Hernandez et al.ll2005h . Right panel: Same as left panel, but in this model we have assumed complete lateral mixing, 
where all of the material within given velocity bins is completely mixed. Two simulations are shown for two different mass cuts, 
Mcut = 2.41 M0 (solid line with open circles) and Mcut = 7.56 Mq (dashed-dotted line with open circles). 



and Ti). The O abundance was derived using the 
O I triplet at 7771-5 A, a robust O abundance in- 
dicator; it appears to be a factor of 4 more abun- 
dant than in the Sun, and more than a factor of 
2 overabundant in comparison with stars of the 
solar neighborhood having similar iron content. 
These photospheric abundances suggest possible 
contamination from nucleosynthetic products in 
the supernova/hypernova explosion that formed 
the compact object in this system. 

The pecuhar velocity of the LMXB V404 Cygni, 
Vpoc = 39.9 ±5.5 kms~^, is larger than the Galac- 
tic velocity dispersion, suggesting that the black 
hole in this system formed in a supernova event. 
Thus, the secondary star could have kept a record 
of this SN/HN event in the chemical composition 
of its atmosphere. We explored this possibility 
using a variety of SN/HN explosion models for 
different geometries. We compared the expected 
abundances in the companion star after contam- 
ination of nucleosynthetic products from initial 
metallicities of the secondary star adopted using 
the average values in solar-type thin-disk stars. 
Metal-rich spherically symmetric models are able 
to reproduce the observed abundances relatively 
well, except for O and Mg, which appear to show 
too high and too low (respectively) abundances in 
comparison to model predictions, regardless of the 



choice of the mass cut. 

A spherical explosion easily provides the energy 
required to explain the peculiar velocity of this 
system, which is in the Galactic plane. There- 
fore, an asymmetric SN/HN explosion may not be 
required. However, we have also modelled the ex- 
pected abundances in the secondary star after cap- 
turing a significant amount of material from an 
aspherical explosion. These models also provide 
good agreement for all elements, as in the spher- 
ical case, except for O and Mg, without invoking 
extensive fallback and mixing. 



In a forthcoming paper (jGonzalez Hernandez et al 



201 iL in preparation) , we will present the Li abun- 
dance and study the ^Li/^Li isotopic ratio to ex- 
plore different evolutionary scenarios and possible 
production of Li due to the X -ray radiation in thi s 
black- hole binary system (see Casares et al.ll2007 ) . 
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